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1  Introduction 


The  Articulated  Total  Body  (ATB)  is  used  at  the  Armstrong  Laboratory  (AL)  for  predicting 
gross  motion  of  the  human  body  under  various  dynamic  environments.  The  model  treats  the 
individual  segments  of  the  human  body  as  rigid  bodies.  However,  the  responses  predicted  by  the 
model  may  not  be  accurate  since  it  ignores  the  effects  of  small  deformations  within  the  individual 
segments  [1-3].  In  the  case  of  human  body  or  dummies,  the  neck  undergoes  large  deformation  which 
may  not  be  modeled  accurately  through  rigid  body  dynamics.  For  example,  ATB  simulation  does  not 
provide  a  good  measurement  of  the  Hybrid  III  Dummy  pendulum  test  [4]. 

The  objective  of  the  paper  is  introduce  the  formulation  of  the  new  ATB  model  which  can 
incorporate  any  number  of  deformable  segments  and  demonstrate  its  advantage  over  the  old  (rigid) 
model  through  simulation  and  comparison  with  Hybrid  III  neck  experimental  test  results. 

It  is  assumed  that  the  deflection  of  a  deformable  body  relative  to  its  own  reference  frame  is 
small  and  linear  such  that  it  can  be  defined  by  a  linear  combination  of  vibration  modes.  In  order  to 
reduce  the  size  and  complexity  of  the  problem,  a  small  number  of  modes  may  be  selected  to 
approximate  the  displacement  field  using  Ritz  approximation  [5].  The  vibration  normal  modes  are 
determined  using  finite  element  modal  analysis  of  the  deformable  bodies. 

First,  the  constrained  equations  of  motion  of  open  chain  multibody  systems  representing  both 
rigid  body  motion  and  small  deformation  are  developed.  Next,  finite  element  modal  analysis  of  the 
Hybrid  III  manikin  neck  is  presented.  The  relevant  frequencies  and  mode  shapes  obtained  from  the 
finite  element  model  of  the  Hybrid  III  dummy  neck  are  used  in  the  ATB  mode!  simulation.  Both 
static  and  dynamic  loadings  are  used  and  the  results  are  compared  with  the  data  obtained  using  the 
Static  Neck  Tester  [6]  and  Head/Neck  Pendulum  [4],  respectively. 
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2  Equations  of  Motion  of  Deformable  Bodies 

To  define  the  gross  motion  of  a  deformable  segment,  a  body  reference  frame,  xyz,  is  attached 
to  the  undeformed  state  of  the  body  and  its  location  (x)  and  orientation  (D)  is  specified  with  respect 
to  an  inertial  reference  frame  AfYZ.  The  elastic  displacement  field  of  a  deformable  body  (a)  may  be 
represented  as  a  linear  combination  of  a  set  of  selected  deformation  modes  determined  by  finite 
element  modeling  and  Ritz  approximation.  Assuming  the  body  is  modeled  with  n  nodes  and  m  modes 
are  selected,  the  displacement  field  containing  nodal  translations  and  rotations  of  all  nodes  is  written 
as: 

a  =  Ta  (1) 


where  a  is  an  mxl  modal  coordinate  vector  and  Y  is  a  6/ixm  modal  matrix.  Therefore  as  shown 
in  figure  1,  the  position  vector  of  a  node  k  on  the  body  may  be  written  as: 


=  X  +  =  X  +  D^(5^+ilr^a) 


(2) 


where  and  are  the  position  vectors  of  nodes  k  in  the  undeformed  and  deformed  states  with 
respect  to  xyz,  is  the  transpose  of  the  direction  cosine  matrix,  and  ijr*  is  a  size  3xm  translational 
modal  submatrix  corresponding  to  node  k.  Taking  the  first  and  second  time  derivatives  of  Eq.  (2), 
velocity  and  acceleration  vectors  of  node  k  are  written  as: 

ijfc  =  X  +  D^(-5j^a)+i|r^d)  (3) 

x^  =  X  +  +  6(cijj^+2\|r^d)]  (4) 
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where  (O  is  the  angular  velocity  vector  about  Jiyz  axes,  and  operating  on  a  vector  produces  a 
skew-symmetric  matrix  for  cross  product  representation  in  matrix  algebra. 

The  variational  equations  of  motion  of  a  deformable  body  at  time  for  a  virtual  displacement 
field  that  is  consistent  with  the  constraints,  are  written  as  [7]: 

-fpdxlx^dV  +  fbxlf^dV  = 

V  V  V 

where  p  is  material  density,  /*  is  the  body  force  at  node  k,  6x*  is  virtual  displacement  of  node  k 
consistent  with  constraints,  V  is  the  undeformed  volume,  and  t*  and  e*  are  defined  as  stress  and  strain 
vectors.  Taking  the  variation  of  Eq.  (2),  the  virtual  displacement  of  node  k  is  written  as: 

=  6x  +  +t|r^6a)  (6) 

where  Sit  is  a  virtual  vector  along  the  axis  of  rotation  and  its  magnitude  is  the  angle  of  rotation. 

For  linearly  elastic  deformation,  the  modal  strain  matrix  relates  strain  and  displacement  and 
the  material  property  matrix  Dg  relates  the  stress  and  strain  vectors  as: 

“C  =  ^0®  =  ^0^^  (7) 

From  Eqs.  (4)  &  (6-7)  and  the  fact  that  the  variational  equations  must  be  valid  for  all 
admissible  6x,  67t,  and  6a,  Eq.  (5)  for  a  system  of  interconnected  bodies  can  be  written  as: 

Af^jc  +  Af^d)  +  (8) 

Mix  H-  +  M^d  >  Blf  ^  Bix  = 


3 


and  the  kinematic  constraint  equations  are: 


Bjjd)  +  +  ^24^  =  v. 


(9) 


Vectors/and  x  are  reaction  forces  and  moments  at  the  kinematic  joints.  The  mass  matrix  elements 
A/^  and  RHS  vectors,  have  been  derived  by  Yoo  and  Haug  [8].  Forces  due  to  strain  energy  and 
modal  damping  forces  are  derived  as: 


4  =  a 

f,  =  d 


(10) 


where  K,  and  Q  are  the  finite  element  stiffiiess  and  damping  matrices  which  are  diagonalized  by 
modal  matrix  Y  in  the  above  equation  [9  ].  The  fiarces  in  Eq.  (10)  are  simply  subtracted  from  the  third 
equation  in  Eq.  (8).  In  the  above  equations,  the  modal  basis  is  reduced  for  the  deformable  segments 
based  on  the  methods  presented  in  Refs.  [10]  and  [11]. 

Submatrices  B~s  are  defined  based  on  the  type  of  kinematic  joints.  The  first  set  of  constraint 
equations  are  defined  for  free  or  ball  and  socket  joints  as: 


X.  -  Xy  -  Djs-Cii  +  DjsjCij  *  = 

-  Dfd).  (6,.  s .  +  2ilr,  d)  +  Z) Jd)y (<bjSj  +  ly^jdp 


(11) 


The  second  set  of  constraint  equations  are  defined  for  pin  joints  as: 


A 


+  D^hh.DjBa-d.  +  Xhh^'Z  = 

J  J  J  rty  «y*^;  ; 

-  Djhj((bjhj(<Dj^2<D„)  +  a)„.A.O)„,  +  A.C„,) 


(12) 


or  for  fixed  joints  as: 


Dj6>.  -  Dj6>i  t  -  DjD,^B,,ii.a.  = 

-  )  *  D[(fi),.u„  4-C.  ) 


(13) 


Matrices  in  the  above  equation  are  described  in  detail  by  Ashrafiuon,  et  al.  [12]. 

The  new  equations  of  motion,  kinematic  constraints,  and  other  required  changes  have  been 
coded  into  the  ATB  model.  The  only  new  input  required  for  the  new  ATB  model  are  the  deformable 
segment  definitions;  e.g,  nodes,  fi'equencies,  etc.  The  new  model  typically  takes  2  to  6  times  the  CPU 
time  of  the  old  model  but  is  Justified  because  of  its  significant  improvement  in  accuracy. 


3  Hybrid  III  Manikin  Neck 

The  Hybrid  III  neck  is  segmented  with  three  3.4  inch  diameter  aluminum  plates  between  Butyl 
rubber  sections  to  simulate  the  human  vertebral  disks,  as  shown  in  figure  2.  The  center  rubber 
sections  are  2.7  inches  in  diameter.  They  are  offeet  0.2  inches  towards  the  front  of  the  neck  and  have 
slits  to  provide  a  different  response  in  flexion  than  in  extension  bending.  There  are  also  aluminum 
end  plates  to  facilitate  assembly  with  a  manikin.  A  steel  cable  runs  through  a  0.625  inch  diameter 


hole  in  the  neck.  The  cable  is  placed  in  tension  by  rotating  an  adjusting  nut  on  a  threaded  portion  of 
the  cable  to  limit  excessive  bending  of  the  neck.  The  total  length  of  the  neck  is  5.66  inches. 

Since  modal  analysis  is  linear,  it  is  not  possible  to  model  the  slits,  the  cable,  or  the  nodding 
blocks  which  are  located  at  the  head/neck  attachment  Joint  To  model  the  slits  in  the  neck  would 
require  using  a  contact  surface  which  implies  a  "force”  as  a  function  of  time.  Since  the  classical 
eigenvalue  problem  does  not  have  a  force  vector,  this  is  not  possible.  The  cable  and  the  nodding 
blocks  produce  nonlinear  effects,  which  again  can  not  be  included  in  the  linear  analysis.  However, 
these  effects  can  be  approximated  by  attaching  a  nonlinear  rotational  spring  at  the  head/neck  joint  and 
using  trial  and  error  approach  to  find  its  parameters. 

The  Hybrid  III  neck  is  bolted  at  one  end  in  the  experimental  tests.  There  are  four  bolts  at  the 
base  and  the  nodes  corresponding  to  these  bolts  are  constrained  in  all  directions.  This  model  can  also 
be  explained  in  terms  of  the  mode  shapes  of  a  cantilever  beam,  figure  3  shows  a  view  of  the  Hybrid 
III  finite  element  mesh;  it  has  2386  elements  and  3427  nodes.  The  material  properties  used  are 
Young’s  modulus,  density,  and  poisson’s  ratio  which  are  respectively  10,000  kpsi,  0.0981  Ib/in^,  and 
0.33  for  aluminum  and  1.2  kpsi,  0.0343  Ib/in^  and  0.49  for  Butyl  rubber. 

The  first  two  modes  (36.1  Hz  and  36.3  Hz)  of  the  Hybrid  III  neck  are  the  first  bending  modes 
of  a  cantilever  beam.  Figure  4  shows  the  first  mode  shape.  The  second  mode  exhibits  some  twist, 
probably  due  to  the  asymmetric  geometry.  The  difference  between  flexion  and  extension  of  the 
physical  neck  is  not  apparent  in  the  finite  element  model;  this  too  has  to  be  adjusted  in  the  ATB 
simulation.  The  third  mode  is  the  first  torsion  mode  excited  at  a  frequency  of  57.2  Hz.  The  second 
set  of  bending  modes  correspond  to  the  fourth  and  fifth  modes  at  firequencies  of  141.0  Hz  and  144.3 
Hz.  Finally  the  sixth  mode  shape  is  the  second  torsion  mode  at  184.3  Hz. 
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4  Static  and  Dynamic  Test  Simulations 

In  the  Static  Neck  Tester  (SNT),  the  Hybrid  III  neck  behaves  similarly  to  a  cantilever  beam 
with  a  load  applied  at  its  free  end,  as  shown  in  figure  5.  The  load  (jF)  is  linearly  increased  from  0  to 
400  lbs  in  2  seconds.  This  simulation  is  "slow  enough"  such  that  dynamic  effects  are  negligible.  See 
Re£  [6]  for  detaUed  description  of  the  experiment.  The  bending  modes  obtained  from  finite  element 
modeling  are  utilized  to  represent  neck  deformation. 

Figure  6  shows  a  comparison  of  the  simulation  results  using  the  new  ATB  model  and  the  SNT 
test  results  reported  by  Spittle,  et  al.  [4]  for  Hybrid  III  static  neck  flexion.  It  can  be  seen  that  the 
bending  stiflfoess  is  predicted  to  be  much  higher  by  the  simulation  for  small  deformation  but  slightly 
lower  for  large  deformation.  This  is  partly  because  dynamic  elastic  properties  of  rubber  are  used  for 
the  model,  which  are  normally  higher  than  the  static  properties,  and  also  nonlinear  (but  still  elastic) 
behavior  of  the  rubber  is  ignored. 

An  ATB  simulation  model  of  the  Head/Neck  Pendulum  (HNP)  test  is  shown  in  figure  7.  The 
HNP  tests  are  dynamic  tests  where  the  head  inertia  may  have  a  significant  effect  on  the  response. 
There  are  three  segments  (pendulum  arm,  neck,  and  head)  and  two  fixed  joints  (jl  &  J2)  connecting 
them.  Therefore,  both  neck  rotation  (4>')  and  head  rotation  ((f))  are  purely  due  to  neck  deformation. 
The  mass  and  geometric  properties  for  the  Hybrid  III  neck  and  head  reported  by  Kaleps,  et  al.  [13] 
are  used.  A  damping  ratio  of  C  =  0.9  is  used  for  the  neck  which  has  an  effective  value  of  about  0.2 
since  the  head  introduces  significant  increase  in  system  inertia.  See  Ref.  [4]  for  a  detailed  description 
of  the  experiment. 

The  excitation  is  provided  through  deceleration  data  obtained  from  the  HNP  tests.  The 
deceleration  data  and  the  initial  impact  velocity  are  applied  along  the  X-axis  of  the  pendulum  arm  for 
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flexion/extension  tests  and  Y-axis  for  lateral  tests.  Several  flexion  tests  have  been  performed  by 
dropping  the  pendulum  arm  from  20“,  40“,  60“,  80“,  and  120“  angles  and  one  lateral  test  by 
dropping  from  65®  resulting  in  impact  velocities  ranging  from  54.74  to  275.4  in/sec.  The  first  four 
bending  modes  which  represent  two  flexion/extension  and  two  lateral  deformation  modes  are 
selected.  These  modes  are  sufficient  for  accurate  modeling  of  neck  deformation  due  to  physical 
characteristics  of  the  system  and  range  of  frequencies  obtained  from  the  deceleration  data. 

The  results  obtained  by  the  deformable  model  were  compared  with  the  rigid  body  model  [6] 
and  the  HNP  test  results.  The  deformable  model  proved  to  be  much  more  accurate  than  the  rigid 
model  consistently  for  all  tests.  This  is  evident  in  figures  8  and  9  where  head  rotation  and  forward 
acceleration  are  plotted  as  obtained  from  the  HNP  20“  test  and  deformable  and  rigid  models.  Note 
the  significant  improvement  in  predicting  the  frequency  of  oscillation  in  addition  to  the  amplitude. 

Figures  11-13  show  comparison  of  head  and  neck  rotation  ((})  and  <j>')  obtained  from  the  new 
ATB  simulations  and  HNP  20“,  60“,  and  120“  tests.  It  can  be  seen  that  the  simulation  follows  the 
test  data  closely  in  most  cases.  However,  for  the  first  two  cases  (20 “  and  60 “)  the  test  data  seems 
to  suggest  that  the  neck  is  oscillating  about  a  negative  mean  value.  This  could  be  due  to  a  defect 
(slight  permanent  defijimation)  in  the  neck.  Also  note  that  the  fundamental  natural  frequency  of  the 
system  is  about  5  Hz  approximately  7  times  smaller  than  neck's  first  (bending)  natural  frequency. 
This,  of  course,  is  because  of  the  addition  of  head  mass  to  the,  neck  which  has  a  similar  effect  on  the 
system  damping  as  observed  by  the  peak  to  peak  ratio.  Finally,  the  sharp  peak  observed  in  the  120“ 
case  (figure  12)  may  be  due  to  the  linear  deformation  assumption. 

Figures  9, 14  &15  compare  the  head  C.G.  forward  (x)  accelerations.  The  simulation  results 
in  all  cases  fijllow  the  same  pattern  as  the  test  data.  However,  the  initial  high  acceleration  peaks  are 
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not  accurately  predicted  by  the  simulation.  These  accelerations  may  be  due  to  activation  of  the 
nodding  blocks  at  the  head-neck  joint  and  the  pre-torqued  cable  whose  effects  are  nonlinear  and  not 
included  in  the  model. 

Figure  16  compares  the  resulting  moments  of  the  60®  flexion  and  65®  lateral  simulation  tests. 
This  comparison  clearly  demonstrates  that  the  high  peak  occurring  at  about  0.05  seconds  into  the 
simulation  is  due  to  the  interaction  of  the  neck  with  the  nodding  blocks  since  the  peak  is  not  observed 
in  the  lateral  test. 

5  Conclusion 

Mathematical  formulation  of  a  new  capability  for  the  ATB  model  has  been  presented  which 
allows  for  modeling  and  integration  of  deformable  segments  into  the  old  rigid  body  formulation.  A 
finite  element  model  of  the  Hybrid  III  manikin  neck  has  been  developed.  Modal  analysis  was 
completed  on  this  model  and  the  first  six  modes  have  been  extracted  for  inclusion  in  the  validation 
of  the  new  ATB  model. 

The  finite  element  model  of  the  Hybrid  III  neck  has  been  incorporated  into  the  ATB  for  quasi¬ 
static  (Static  Neck  Tester)  and  dynamic  (Head/Neck  Pendulum  test)  simulations  and  comparison  with 
experimental  data.  The  model  was  shown  to  be  stiffer  than  the  actual  neck  in  the  static  simulations 
because  only  dynamic  elastic  properties  were  used.  Dynamic  simulations  showed  the  deformable 
model  to  be  much  more  accurate  than  the  rigid  model. 

Overall  the  new  deformable  element  modeling  methodology  works  well  in  the  ATB  model  and 
can  be  applied  to  model  other  deformable  structures  such  as  a  seat  back  or  the  human  neck  and  spine. 
The  Hybrid  III  neck  and  other  deformable  models  will  be  used  in  whole  body  simulations  to  fully  test 
and  validate  the  new  capability. 
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Cable 


Figure  2.  Hybrid  EH  manikin  neck 
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Figure  3.  Hybrid  III  Hnite  element  mesh 
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Figure  4.  Hybrid  III  first  mode  shape  (36.1  Hz) 
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Figure  7.  Head/Neck  Pendulum  test  model 
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Figure  9.  Comparison  of  Hybrid  IH  head  rotations  in  B[NP  20°  flexion  test 
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Figure  10.  Hybrid  HI  rotations  in  HNP  20®  flexion  test 
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Figure  12.  Hybrid  ffl  rotations  in  HNP  120®  flexion  test 
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Figure  13.  Hybrid  HI  rotations  in  HNP  65*  lateral  test 
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Figure  14.  Head  x  acceleration  in  HNP  60®  flexion  test 
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Figure  15.  Head  x  acceleration  in  HNP  120**  flexion  test 
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Figure  16.  Moments  for  60**  flexion  and  65**  lateral  tests 
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